Dynamics of the entanglement between two oscillators in the same environment 
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We provide a complete characterization of the evolution of entanglement between two oscillators 
coupled to a common environment. For initial Gaussian states we identify three phases with different 
qualitative long time behavior: There is a phase where entanglement undergoes a sudden death (SD) . 
Another phase (SDR) is characterized by an infinite sequence of events of sudden death and revival 
of entanglement. In the third phase (NSD) there is no sudden death of entanglement, which persist 
for long time. The phase diagram is described and analytic expressions for the boundary between 
phases are obtained. Numerical simulations show the accuracy of the analytic expressions. These 
results are applicable to a large variety of non-Mar kovian environments. The case of non-resonant 
oscillators is also numerically investigated. 



Entanglement is responsible for the most counterin- 
tuitive aspects of quantum mechanics. In the early 
days of quantum physics, entanglement motivated many 
philosofi-cal discussions but nowadays it is viewed merely 
as a source of surprises but, mainly, as a physical re- 
source. In fact, entanglement is the key ingredient for 
quantum teleportation and is believed to be the origin 
of the power of quantum computers. Systems of qubits 
interacting with common or independent environment 
where extensively analyzed in the literature In 
such context, the existence of a peculiar property was no- 
ticed in Entanglement, as opposed to energy, can dis- 
appear from a system in a finite time (this phenomenon 
was denoted as "sudden death" SD of entanglement). 

Here we study the evolution of the entanglement be- 
tween two oscillators interacting with a common bosonic 
environment. Entanglement dynamics for such kind of 
systems was analyzed before and a variety of results 
is available in the literature. For example: It was 
shown that the environment may completely dissentan- 
glc initially entangled state (a two-mode squeezed state). 
Thus, in 0, 0] the existence of sudden death was 
proved and analyzed for such case. The evolution of the 
same type of initially entangled state was analyzed when 
the two oscillators interact with a common bath under 
a simplifying assumption: the Markovian approximation 
[Sj JH, JjJ- In that case, a condition for the existence of 
an asymptotic entangled state was deduced 0]. If such 
condition is not satisfied (see below) sudden death takes 
place. Also, it was shown that the interaction with a 
common environment opens the door to the creation of 
an entangled state from an initial separable state 



Most of these works were done under some sort of Marko- 
vian approximation. More recently, the non-Markovian 
regime was analyzed 



12, 13, [li|. Under special condi- 



tions a simple result emerged: The final state remains 
entangled but half of the initial entanglement is lost [3] . 

Our work provides a unified picture to understand the 
origin of the different qualitative behaviors ("phases") 
that characterize the evolution of entanglement for long 
times. We will fully characterize these phases of the evo- 



lution and obtain equations describing the boundary be- 
tween them. For this purpose we will use a well known 
tool: the exact master equation for quantum Brownian 
motion [l^. We consider two harmonic oscillators {xi 
and X2) bilinear ly coupled between them and with a col- 
lection of harmonic oscillators 15, 3, 12, The total 



Hamiltonian is H — Hs + Hint + H^nv where 
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Using coordinates x± = {xi ± X2)/V^, the Hamilto- 
nian transforms into Hs — {p% + p^)/2m + ra{L0^x'^_ + 

'^+a;^)/2 -I- TOC-i x+x^, where w^. — [uol -f a;|)/2 ± C12 

and c-i = {ujI — w|)/2. This model can be exactly 

solved assuming that the initial state of the environment 
is thermal with initial temperature T [l5| . and entangle- 
ment properties between the system and the environment 
were studied in The resonant case (wi = a'2), where 
the x± oscillators are decoupled, is the simplest: The ex- 
act master equation for the reduced density matrix of the 
two oscillators p is (S. = 1) [H, [3]: 



1 



[Hfi,,p] - i"f{t)[x+,{p+,p}] 



D[t)[x+, [x+,p]] - f[t)[x+, [P+.P]_ 



(2) 



Here, the renormalized Hamiltonian is Hn ~ Hs + 
m5uj'^{t)x\/2. The coefficients 5uj'^{t), -/{t), D{t) and 
f{t) depend on the spectral density of the environment. 



defined as J(w) — X]„ c^<5(a; 



J/2m„ 



and also on 



its initial temperature. The explicit form of these coef- 
ficients is rather cumbersome and was studied in detail 
elsewhere 15, 2^. We will focus on the ohmic environ- 
ment with a spectral density J{u}) ~ 2m7ow0(w — A)/7r. 
The high frequency cutoff A defines a characteristic 



2 



timescale A"-'^ over which the coefficients j{t) and 5uj'^{t) 
vary. For times t ^ A""'^ these two temperature indepen- 
dent coefficients settle into asymptotic values: j{t) 
7 = 270 and Su!'^{t) — > —AmA'^/Tr. The time dependent 



frequencies fll — ' 



Suj'^{t)/2 approach cutoff in- 



dependent values only if the bare frequencies have an 
appropriate dependence on the cutoff. The coupling con- 
stant ci_2 must also be renormalized in the same way so 
that the time dependent coupling Ci2{t) = ci2 + Suj'^{t) /2 
approaches a finite cutoff independent value. The be- 
havior of the diffusion coefficients D{t) and f{t) is more 
complicated and depend on the initial temperature. For 
the moment we just need to mention here that for realis- 
tic environments these coefficients approach asymptotic 
values after a temperature-dependent time (which equals 
A~^ only in the high temperature regime). A word on 
notation: upper case letters will be used to denote renor- 
malized quantities. The time label will be omitted when 
referring to the asymptotic value of the corresponding 
function (i.e., fli^2 denotes the asymptotic value of the 
renormalized frequency of the oscillators, etc). 

From the master equation ^ we can derive simple 
evolution equations for the covariance matrix Vij (t) = 
Tr{p{t){r,,rj})/2 - Tr{p{t)ri)Tr{p{t)r-j) where i,j = 
1, . . . ,4 and r = (a;_,p_,a;+,p+). Some of these equa- 
tions are particularly illuminating: In fact, equations for 
the covariances split into two blocks of 2 x 2. The evolu- 
tion of the first block formed with the second moments of 
X- and p_ corresponds to a free oscillator with frequency 
ri^(t). The evolution of the second block, formed with 
the second moments of x+ and satisfys: 
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These equations contain most of the necessary informa- 
tion to fully analyze the evolution of the entanglement 
between initial Gaussian states. To solve them exactly 
we need to know the time dependent coefficients that ap- 
pear in the master equation. But, remarkably enough, 
we can use the above equations to understand the qual- 
itative behavior of entanglement. For this, we only need 
to assume that the time dependent coefficients in ^ ap- 
proach asymptotic values. In such case, there is a stable 
stationary solution where the dispersions A^x+ = (x^) 
and A^p+ = (p^) are 



Ap+ = 



nAx+ = 



D 
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Entanglement for Gaussian states is entirely deter- 
mined by the properties of the covariance matrix Vij. 
In fact, a good measure of entanglement for such states 



is the logarithmic negativity £V 2l|, I24I , computable as 



2l|,[2l|23| 



E^r = max{0, - ln(2j/,nin)}, 



(6) 



where Vniin is the smallest symplectic eigenvalue of the 
partially transposed covariance matrix. There are known 
expressions for Ej^ for particularly relevant Gaussian 
states which will be used as initial conditions in our study. 
For this reason it is useful to mention them here: For the 
two-mode squeezed state, obtained from the vacuum by 
acting with the operator exp(— r(a|a2 — aia2)), we have 
£V = 2\r\. For this state the dispersions satisfy the min- 
imum uncertainty condition Sx-^-Sp+ = Sx-Sp- = 1/2. 
The squeezing factor determines the ratio between vari- 
ances since m^l6x+/Sp+ = dp-/{mil,6x-) = exp(2r). As 
r ^ 00 the state becomes localized in the_p+ and 
variables approaching an ideal EPR state [2J|. Another 
initial state we will consider is a separable squeezed state 
for which mnSxi^2/Spi,2 — exp(2r). 

The evolution of entanglement for initial Gaussian 
states can then be simply analyzed by using the previ- 
ous results. Let us consider a general initial state. For 
the resonant case, the x- virtual oscillator decouples ex- 
actly. Using the equations Q we see that Aa;+ and Ap+ 
approach asymptotic values after a timescale I/7, fixed 
by the dissipation rate. After that time, the covariance 
matrix in the (a;+,a;_) bases has a 2 x 2 block with os- 
cillatory functions (corresponding to the X- oscillator) 
and another 2x2 block, corresponding to the x+ vir- 
tual oscillator, which is diagonal. Using this form for 
the covariance matrix (and changing basis to obtain co- 
variances of the original xi^2 oscillators) it is simple to 
find its smallest symplectic eigenvalue and compute the 
following logarithmic negativity [2^ 



Ej^{t) max{0, E^ + AE^G{t)}. 



(7) 



Where G{t) is a function with period ir/il^ in { — 1, 1}. 
The mean value Ej^f about which Ej^ oscillates is 

Em ^ max{r, Vcrit} - ^ ln(2Aa;+ Ap+), (8) 



and the amplitude of the oscillation is: 
AE^r = min{r, Vcrit}- 



(9) 



Here r is the initial squeezing factor defined as 
r = |i ln(m57_|^)|, rcrit is related to the squeezing fac- 
tor of the equilibrium state for the x+-oscillator defined 
as rcrit = I5 lii(™^-s^)l- The dispersions Ax+ and 
Ap^ are the asymptotic values ([5]). 

These simple results will enable us to draw gen- 
eral conclusions about the dynamics of entanglement 
for long times. Three qualitatively different scenar- 
ios emerge: First, entanglement may persist for arbi- 
trary long times. This phase, which we call "NSD" 
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FIG. 1: Phase diagram for ohmic environment (il = 1, 
70 = 0.15, A = 20, m = 1, C12 = 0). The sudden death 
(SD), no-sudden death (NSD) and sudden death and revival 
(SDR) phases describe the three different qualitative long 
time behaviors for the entanglement between two oscillators 
interacting with the same environment. The SDR phase is 
centered about the dashed line Sr and has a width given by 
the dotted line rcrit- This is the case for temperatures above 
the one for which 5*,. — rcrit- Below this temperature the 
role of Sr and rcrit are interchanged. SDR separates the SD 
and NSD phases. The low temperature NSD island is due to 
non-Markovian and non-perturbative effects. Ej^ in the NSD 
phase is the distance to the dashed line for r > rcrit, and the 
distance between the dashed and dotted lines for r < rcrit- 



(for no-sudden death), is realized when the initial state 
is such that £V — ^Ej\f > 0, which translates into 
|f — rcrit\ > ln(2Aa;+Ap+)/2 = Sr- Then, there is 
a phase where entanglement undergoes an infinite se- 
quence of events of "sudden death" and "sudden revival" 
1261. I27II. This occurs if the initial state is such that 



|ln(l/2ml]_A2a:+)/2| 



< 



< 



ln(2A2p+/mrj_)/2. We 



denote this phase as "SDR" (for sudden death and re- 
vival). Finally, a third phase characterized by a final 
event of "sudden death" of entanglement may be real- 
ized if r < ~Ec = ~ln(l/2TOrj_A2x+)/2. This phase is 
simply denoted as "SD" (for sudden death). 

Depending on the properties of the environment (ini- 
tial temperature, damping rate, etc) a given initial state 
(parameterized by the squeezing r) will belong to one of 
the three phases. For the ohmic environment we can use 
exact analytic expressions for the coefficients D and / 
(20| to obtain the phase diagram displayed in Fig. [TJ In 
the diagram, the areas corresponding to each of the three 
phases are displayed. As a reference, we also include the 
curves that show the temperature dependence of Sr and 
rcrit (dashed and dotted lines respectively). 

The phase diagram provides complete information 
about the asymptotic behavior of entanglement. It is 
worth mentioning some features of the diagram. The 
NSD phase present at low temperatures is purely non- 
Markovian and non-perturbative. Its area shrinks as the 
damping rate decreases. This phase includes the coher- 
ent states, that can become entangled at very low tem- 
perature. The zero temperature line is also interesting: 
It contains states in the NSD phase for small and large 
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FIG. 2: Logarithmic negativity for resonant oscillators in the 
same environment, (a) For T = the NSD phase appears 
both for for large and small squeezing. Asymptotic behavior 
of initially entangled or separable states only depends on r. 
The amplitude of oscillations vanishes when r — > 0. (b) The 
SDR phase appears for intermediate values of squeezing at 
zero temperature, (c) T/Q, = 10, the SD phase appears for 
small r and NSD phase for large squeezings, oscillations in 
the steady state are attenuated as the temperature increases. 



squeezings. However for an intermediate range of squeez- 
ings centered about rent states belong to the SDR phase. 

To the contrary, the high temperature region of the 
diagram is quite different. Thus, for high temperatures 
we have Ec < (which implies that coherent states do 
not get entangled) and also rcrit ~^ (which implies 
that the region covered by the SDR phase becomes nar- 
rower). Thus, for high temperatures initial states with 
large squeezing (r > ln(2Aa::+Ap+)/2 = Sr) retain some 
of their entanglement while those with squeezing fac- 
tors smaller than the critical value Sr suffer from sudden 
death. We remark that this is almost the same condi- 
tion obtained previously for persistence of entanglement 
in the Markovian regime ^i^]. However, our analysis shows 
that the boundary between SD and NSD phases is rather 
subtle: for any finite temperature the two phases are sep- 
arated by a very narrow portion of SDR phase (in this 
phase there are oscillations of the entanglement whose 
amplitude, rcrit decreases as temperature grows). This is 
yet another interesting non-Markovian effect identified by 
our analysis. A final comment on the phase diagram: The 
NSD phase is characterized by a non-vanishing asymp- 
totic entanglement that can be quantified in a straight- 
forward way from the phase diagram itself. The average 
value of the logarithmic negativity is simply the distance 
to the dashed line (which signals the midpoint of the 
SDR phase) or just the distance between the dashed and 
dotted lines for r < rcrit - 

The validity of the above rather simple but nontrivial 



analytical predictions can be verified by an exact numer- 
ical solution to the problem. We considered an ohmic 
environment with the following parameters: 70 = 0.15, 
A — 20, m — 1, U — 1, C12 = (extension to the 
case where the oscillators interact can be easily done). 
The results for the evolution of are shown in Fig. 
We clearly see that the final entanglement produced by 
different initial states only depends upon the squeezing 
factor r. The existence of events of sudden death (and 
revivals) can also be seen from the numerical solution, 
and is predicted from our analytic results characterizing 
the SDR phase. Fig. [21 It is also noticeable that the am- 
plitude of the oscillations vanishes when //D ^ 0, which 
is a feature of the high temperature limit. The oscillat- 
ing behavior of the final entanglement is a consequence 
of the non-zero asymptotic value of the anomalous diffu- 
sion coefficient f{t), (see Fig. [2]). Moreover, this same 
coefficient is also responsible for the generation of asymp- 
totic entanglement in initial coherent states. But in this 
case the entanglement is constant Ej\f = Ec and increases 
with 7o, as the final state of the oscillator x+ becomes 
more and more localized in position. This is another 
purely non-Markovian effect that is entirely lost in the 
RWA-Markovian approximation With the numerical 
solution, we verified that every possible temporal evolu- 
tion that can be obtained by changing parameters fits 
either in the NSD, SDR or the SD phase. The bound- 
aries between phases is also accurately described by the 
conditions mentioned above. 

The above properties are valid under a single impor- 
tant assumption: the two oscillators are resonant. If this 
is not the case the analysis becomes more complicated: 
the master equation is no longer valid since there are 
terms coupling the a;_ and x+ virtual oscillators. Due to 
this coupling the a;_ oscillator is not protected from the 
environment and also reaches an equilibrium state. In 
Fig. [3] we show how things change when we move away 
from the resonance condition. It is clear that Ej\f decays 
much faster in this case. Also, the value of E^ was ana- 
lyzed away from resonance for different times. The result 
is shown in Fig. [3] where it is clear that the resonance 
peak becomes sharper and sharper as time grows. 

Our results provide a complete picture of the possible 
behavior of entanglement for resonant oscillators coupled 
to a common environment. The same analysis is valid for 
any environment such that an asymptotic value of the 
coefficients of the master equation is obtained. For ex- 
ample, for super-ohmic environments the damping coef- 
ficient approaches a small value while both f{t) and D{t) 
vanish asymptotically. Our analysis predicts that there 
will not be an equilibrium value for the dispersions of the 
x+ oscillator and that entanglement will oscillate approx- 
imately recovering its initial value after some time. This 
is indeed observed in our numerical simulations (contrary 
to what is indicated in 13|, where the results obtained 
seem to be an artifact produced by an unreasonably low 
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FIG. 3: Entanglement dynamics for non-resonant oscillators 
initially in a squeezed separable state (r = 3) and T /Q. = 10. 
(a) Entanglement is created between non-resonant oscillators 
but it banishes in finite time. (&) Ej^ for different times as a 
function of the frequency of the second oscillator, the resonant 
condition is essential for asymptotic entanglement. 



cutoff). Finally, we remark once again that our results for 
the asymptotic entanglement are highly dependent on the 
validity of the resonance condition. For non-resonant os- 
cillators the generic fate for long times is disentanglement 
due to the interaction with the environment. This issue 
will be analyzed in detail elsewhere. Authors acknowl- 
edge support from CONICET and Anpcyt (Argentina). 
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